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ABSTRACT 



We provide a complete theory of the phase closure of a binary system in which a small, feeble, and unresolved companion acts as a 
perturbing parameter on the spatial frequency spectrum of a dominant, bright, resolved source. We demonstrate that the influence of 
the companion can be measured with precision by measuring the phase closure of the system near the nulls of the primary visibility 
function. In these regions of phase closure nulling, frequency intervals always exist where the phase closure signature of the compan- 
ion is larger than any systematic error and can then be measured. We show that this technique allows retrieval of many astrophysically 
relevant properties of faint and close companions such as flux, position, and in favorable cases, spectrum. We conclude by a rapid 
study of the potentialities of phase closure nulling observations with current interferometers and explore the requirements for a new 
type of dedicated instrument. 
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interferometers - methods: observational - techniques: interferometric - astrometry 



1. Introduction 

A common feature of the observational problem posed by spec- 
troscopic binaries, evolved stars companions, hot Jupiters, and 
extrasolar planets is the presence of a point-like source of low 
flux in the vicinity of a much larger and much brighter star. 
1 Indeed, it is the properties of the faint companion that are de- 
' sired, and are the most difficult to acquire. Directly detecting the 
' flux of an extrasolar planet and establishing its spectrum seems 
the most demanding task, to be fulfilled only by dedicated spa- 
tial missions (DARWIN: Cockell et al. 2008, and TPF -I: Lawson 
6t ^1- 2008) using interferometry nulling techniques dBracewelll 
i 119781) . 

Ground-based optical long baseline interferometers also 
have the capability of detecting extrasolar planets restricted to 
' the case of a gian t planet orbiting close to its parent stars 
dVannier et alj|200^ . such as 5 1 Peg th e first ex oplanet to be un- 
raveled by radial velocity techniques dMavor & Oueloz 1995t). 
The principle consists in measuring interferograms in different 
spectral regions where the ratio between the planet photons and 
the stellar ones changes because of an emission line or an ab- 
sorption line in the planet atmosphere. The position of the inter- 
ferogram measured by the phase is directly related to the posi- 
tion of the photocenter of the star-planet system on the sky. The 
slight shift in the phase between two spectral regions is in princi- 
ple measurable for bright extrasolar planets such as hot Jupiters. 

Even when using color-differential techniques the phase 
measurement is plagued by the variation of the refractive index 
in air. This is why efforts are currently concentrated on using 
the phase closure, which corresponds to the sum of the phases 
measured in three connected baselines (see Sect. |2.4| for details). 
This technique allows the observer to remove any atmospheric 
perturbation or instrumental-based errors localized in each arm 



of the interferometer. This detection has not been performed yet 
because of instrumental biases and the need to decrease the pho- 
ton noise even with large telescopes. 



An important property of stars resolved by long baseline in- 
terferometry is that the coherence of the light decreases with 
increasing spatial frequencies down to zero before increasing 
again following the well kn own behavior of the Bessel functions. 
iMichelson & Pease (1921) used this property to measure the di- 
ameter of Betelgeuse for the first time. Indeed, the visibility of 
a photosphere of uniform brightness distribution over a circular 
disk is canceled out at the spatial frequency 0.61//?* where 
is the stellar angular radius. Another property is that the phase 
of the visibility jumps by 180 degrees at each crossing of the 
Bessel function's zeros (hereafter called visibility nulls or sim- 
ply nulls). At each null, the coherent flux of the star is canceled 
out, making it easier to detect a change in the phase jump due to 
a faint close companion. 



In this paper we demonstrate that, in precisely the restricted 
case of the lobe crossing, the influence of even a very faint com- 
panion can be measured with precision, and its spectrum can be 
obtained by measuring the phase closure of the system near the 
nulls of the primary visibility function. In these regions where 
the phase closure nulling of the primary is effective, the phase 
closure signature of the secondary stands out to a point where its 
basic properties (flux, position, spectrum, etc.) can be measured. 
In a companion paper (Duvert etal. 2009), this technique was 
used to detect the companion of the single-lined spectroscopic 
binary crPuppis, at a few stellar radii distance. 
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Fig. 1. Visibility for a star well represented by a uniform brightness distribution over a disk of angular radius /?+(left panel), as a 
function of the reduced frequency (defined as the frequency times the stellar radius). The visibility is often presented under the form 
of amphtude (center panel) and phase (right panel). 



2. Description of tlie problem 

2.1. Visibility of a stellar photosphere of uniform brightness 

When observing a star with a multi-aperture interferometer, it 
is well known that one detects the amplitude of the visibility 
function (or the amplitude squared) and that this amplitude de- 
creases with increasing spatial frequency u - B/A, where B is 
the baseline and A the wavelength. Thanks to the Zernike-Van 
Cittert theorem, it can be shown that this visibility function co- 
incides with the Fourier transform of the spatial energy distri- 
bution of the star; hence, for a star whose energy distribution 
is well-represented by a uniform disk, the visibility function is 
proportional to a Bessel function of the first order: 



VJu) = 2 



J i{2n u Ri,) 



(1) 



The visibility function is a complex number that, in this special 
case of a simple centro-symmetric brightness distribution, turns 
out to be real (i.e. the imaginary part is zero). It is often displayed 
by its amplitude and its phase (see Fig.[Tll. 

The visibility becomes negative when the baseline is longer 
than 0.61 A/Ri,, corresponding to the first null, so the phase 
suddenly changes from zero to 180 degrees (see right panel of 
Fig. HI)- It becomes positive again when crossing the second null, 
and so on. 



2.2. Visibility of a resolved star with a faint companion 

When an object is composed of two components, the object vis- 
ibility function is the sum of the visibility function of the two 
components weighted by the fraction of flux of each component 
and a phase modulated by the apparent angular separation on the 
sky. If one adds a faint companion to the star of the previous sec- 
tion with a flux ratio r and separated from the star by s (assuming 
for sake of simplicity that the system orientation is parallel to the 
frequency axis), then the complex visibility of the system is 



i(u) = 



V^ju) + re- 
1 +r 



i In us 



(2) 



If the flux ratio r is small, the amplitude of the visibility func- 
tion will be very slightly modified by the presence of the faint 
companion. The stronger effect occurs around visibility nulls of 
the primary, where the system visibility is no longer zero. At 



these locations, the visibility perturbation is equal to the flux 
ratio. This effect still remains weak and is beyond the perfor- 
mances of current interferometers as soon as the flux ratio is 
less than about 1%. However, besides the fact the companion 
is not detectable, studies around visibility nulls keep their in- 
terest for precise diame ter and limb darkening measurements 
(iLe Bouquin et al.ll2008l) . 

2.3. Phase changed by the faint companion 

More interesting is the phase of the object visibility. By defini- 
tion, the tangent of the phase of a complex number is equal to 
the ratio of its imaginary part by its real part. Therefore we have 



tan (p{u) - 



r sin(27rMs) 



Vi,(u) + r cos(2nus) 



(3) 



For small flux ratios and except for around nulls, the phase of 
the object is very close to or 180 degrees, and may be approx- 
imated by 

^ ^ rsin(2;rMs) 
0(m) * , ' , + nn (4) 
VJu) 

with n = if Vi,{u)>0, and n = 1 if Vi,{u)<0. For separations 
s/R^, much greater than 1, it is quasi periodic, with a pseudo 
period s 

The phase is the result of two contributions: from the star 
Vi,{u) and from the companion re'^'"'\ The contribution of the 
star decreases as the frequency increases and is even canceled at 
the nulls, while that of the companion keeps constant as long 
as it is not resolved. Most of the time, the companion pro- 
duces a phase signature in the range ±r/y*(M). But more impor- 
tant, around visibility nulls, in the frequency ranges for which 
|V*(m)| < r, this phase signature becomes significant, with an 
exact value of Inu^s at the frequencies mq of the nulls {xns/Ri, 
for the first null), much greater than 180 degrees. It follows that, 
as opposed to the visibility amplitude, even for small flux ratios, 
frequency intervals always exist around nulls, where the phase 
signature of the companion is larger than any systematic error 
and can then be measured. 



2.4. Phase closure of a resolved star with a faint companion 

An interferometer formed by 3 telescopes transmits 3 frequency 
systems, ui2, M23, and uu - un + M23, associated with the 3 
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Fig. 2. Phase closure modulo In of a double system, formed 
by an extended uniform disk and a point source, as a function 
of the maximum reduced frequency. The spectral resolution is 
K = 1500. The 3 frequencies {m23, uii, wn) (see Eq.|5]l have been 
chosen in the ratio 1,2,3, with the maximum frequency around 
the first zero of visibility of the primary star. The phase closure 
is displayed for 3 separations (10/?*, 100/?*, and 1000/?*) and 
a flux ratio of 0.01. The signature of the secondary source is 
dominant in the region around the minimum visibility of the pri- 
mary (near reduced frequency 0.61). The thin line with a 180 
degree shift at 0.61 corresponds to the phase closure of the pri- 
mary alone. 



baselines, each characterized by a phase and a visibility. It pro- 
vides one phase closure defined as the phase of the bispectrum 
/(mi2, M23), with 



/(mi2,M23) =< i(mi2)i(M23)i'(mi3) > 



(5) 



where * denotes the complex conjugate, and <> represents an 
ensemble average. The phase closure has two important proper- 
ties: it is insensitive to atmospheric phase fluctuations, and for a 
point source, it is equal to zero. Then, for an extended object, it 
is given by 

Cf)c = 0o(mi2) + <Po(U23) - <Po{un) (6) 

where (po is the phase of the object spatial Fourier transform. 

If the telescopes are aligned with baselines parallel to the 
direction of the double system, the phase closure may be written 

rsin{2nu\2s) rsin{2nu2T,s) rsin{2nu\is) 
(pc ~ 1 ^ — + 

V*(mi2) V*(M23) V*(mi3) 

("12 + "23 -n 13)^- (7) 

Figure|2]shows the expected phase closure, around the first vis- 
ibility null, from a double system with a flux ratio of 1% and 
various separations. The system is observed with a spectral res- 
olution of 1500, li ke the one available with the AMBER instru- 
ment at the VLTI (iPetrov et al.ll2007l) . Note the importance of 
the phase closure signature from the companion. 
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Fig. 3. Absolute value of the useful phase closure signal of a 
double system formed by a uniform disk and a point source as a 
function of the highest frequency located around the first visibil- 
ity null of the primary. The spectral resolution isK - 1500. The 
phase closure of the uniform disk has been subtracted for clarity. 
From top to bottom, the flux ratio varies from 10"^ to 10"^ and 
3 separations are considered: 10/?* (full lines), 100/?* (dashed 
lines), 1000/?* (dotted lines). The horizontal lines are set to a 
given systematic error of 0. 1 degree, showing the regions where 
the phase closure signature of the secondary becomes measur- 
able. 



within which the useful signal is stronger than any systematic 
error, hence can be measured. 

Figure [3] shows the useful phase closure signal for various 
flux ratios and distances of 10, 100, and 1000 stellar radii. The 
horizontal lines correspond to a given systematic error of 0.1 
degree. The frequency interval within which the useful signal is 
stronger than this error diminishes rapidly with the flux ratio, but 
is always present. 

We propose using the regions around minima of visibility, 
where the phase closure nulling of the primary is effective, and 
the useful signal is over the limit crj,,j imposed by systematic 
errors, to detect and to characterize faint companions. By doing 
so, we do not cancel the stellar flux as in classical nulling exper- 
iments dBracewelll 1 9781) . but we only cancel the stellar coherent 
flux. The main limiting noise will then be the photon noise from 
the central star, as shown in Sect. 13.21 



3. Phase closure nulling 

In the case of a system formed by an extended uniform disk and 
a point source, the useful phase closure signal is the departure 
from or 180 degrees. We have then clearly interest in selecting 
the working frequency interval at the highest frequencies, around 
visibility nulls at m = mq, where the phase closure signature of the 
companion is at its maximum. In these regions, there is always 
(as the stellar visibility crosses the zero) a frequency interval 



3.1. Spectroscopy of close companions 

The spectroscopy of close companions may be extracted from 
spectrally dispersed phase closure measurements obtained with 
an A^tei telescopes interferometer (A^tei^3). The baselines should 
be such that the longest one corresponds to a null of the primary 
star in order to maximize the useful signal. As an example, an 
interferometer with 3 telescopes and frequencies {M23, M12, mo) in 
the ratio 1:2:3, allowing bootstrapping, would be a very efficient 
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system. The data consists of spectrally dispersed phase closure 
measurements performed as a function of the spatial frequency. 

The double system is characterized by 4 parameters: the stel- 
lar radius the flux ratio r, the separation s, and the position 
angle. The 4 parameters may be extracted from a modeling of 
the phase closure variations as a function of both the spatial fre- 
quency and the wavelength. The spectrum of the companion is 
then obtained by multiplying the flux ratio by the spectrum of 
the primary. 

3.2. Error analysis 

We consider for simplicity an interferometer formed by A^tei 
aligned telescopes in a non-redundant configuration. We also as- 
sume that the direction of observation varies little during the ex- 
periment so that the double system is characterized by only 3 
parameters: the stellar angular radius 7?*, the flux ratio r, and 
the projected separation s along the direction of observation. 
Spectrally dispersed phase closure measurements are performed 
as a function of the spatial frequency and then fitted, with a 3- 
parameter model. 

The limiting noises in optical interferometry from the vis- 
ible to the near infrared are the atmospheric noise, the photon 
noise from the source, and the detector readout noise. General 
expressions of the errors on the parameters of the double system 
are derived in Appendix |A] In the photon noise regime, the er- 
rors reach their minimum values and then stay constant as soon 
as the working frequency domains cover at least one period 
For separations greater than a few stellar radii and working fre- 
quency intervals over i the error on the flux ratio and that on 
the separation, with a 3-telescope interferometer, are given by 



air). 
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1 / 0.61 \ 



(8) 
(9) 



where K is the total number of photoevents detected during the 
whole observation and M,„a,i the average frequency of the longest 
baseline. 

These formulae have been checked numerically through a 
rigorous error analysis. They are valid everywhere, even around 
the nulls. Besides the error on the separation being inversely pro- 
portional to the maximum frequency, the errors are independent 
of the working frequency domains. This clearly reinforces the 
usefulness of observing around visibility nulls, where the useful 
phase closure signal is dominant. The error on the flux ratio only 
depends on the total number of collected photoevents, while that 
on the distance is also inversely proportional to the flux ratio. In 
the photon noise regime, both errors are independent of the spec- 
tral resolution "R, as long as the spectral averaging is performed 
over less than half a period 1/2 s, that is, when Umaxfi^^ <0.5 i"' . 
This condition provides a maximum recoverable distance im- 
posed by the spectral resolution of 



-Rxi 



(10) 



The error on the flux of the companion is 
o-(rK) ^ 3 VZ. 



The direct photometric detection of the companion would pro- 
vide an error of yfrK. In terms of performances, the direct detec- 
tion is3/ ^/r better than detection from phase closure. However, 
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Fig. 4. a (top): Error on the flux ratio of a double system 
from the closure phase information of a 3-telescope interferom- 
eter. The signal is made by 900 phase closure regularly sam- 
pled within the working frequency intervals, and spectrally dis- 
persed along 270 spectral channels, with a spectral resolution of 
1500. The experimental parameters are; observing wavelength: 
1.65;um; optical bandwidth: 0.3;um; total integration time: 3 
hours; integration time per frame: 0.2s; detector readout noise: 
lOe ; 32 pixels per frame. The calculations took the photon noise 
and the detector readout noise into account. Three cases were 
considered: 3 telescopes with diameter D=2m, a Strehl ratio 
S=0.5, and 1% total transmission (upper curve, this case cor- 
responds roughly to the present state of the AMBER experiment 
on the VLTI), D=2m, S=0.9 and 10% total transmission (middle 
curve), D=8m, S=0.5 and 10% transmission (lower curve). The 
broken lines correspond to an integration time per frame of 12 
seconds. Represented a re: Siriu s with a possible M dwarf com- 
panion (Benest & Duve njl995 ), the giant K 5 star cr Puppis with 
its G2 dwarf companion (iDuvert et al.ll2009h . and two stars with 
known planetary companions, b (bottom): Baseline needed to re- 
solve the primary star in the H band at 1 .65fim, namely to reach 
the first visibility null, for late dwarfs and giants. The horizontal 
line at 200 m corresponds to the maximum baseUne of the VLTI. 



this is the price to pay because, to our knowledge, there is no 
other direct method than the one we propose, which is capable 
of detecting stellar companions within an Airy disk and at a dis- 
tance from the hosting star as small as a few stellar radii. Once 
the conditions above are fulfilled, measuring the flux ratio for 
each spectral element at resolution K is an eff'ective way to mea- 
sure the companion spectrum. 



(11) 3.3. Performances 



The mean error on the flux ratio in the H band at 1.65//m from 
a 3-telescope interferometer and an integration time of 3 hours, 
obtained by averaging all the spectral channels over an optical 



A. Chelli, G. Duvert, F. Malbet and P. Kern: Phase closure nulling 



5 



bandwidth of 0.3//m, is plotted in Fig.|4^. The upper curves cor- 
respond to 3 telescopes of 2 m, 1% transmission, and a Strehl ra- 
tio of 0.5. It roughly represents the present state of the AMBER 
instrument on the VLTI with Auxiliary telescopes. The middle 
and lower curves correspond to telescopes of 2 m and 8 m, a 
transmission of 10%, and Strehl ratios of 0.9 and 0.5, respec- 
tively (see the caption of Fig.|4]for details). 

If we assume a signal-to-noise ratio of 10 on the flux ra- 
tio for a positive detection, then the AMBER instrument with 
Auxiliary telescopes would allow detection of companions 10"* 
fainter around second-magnitude stars and 10^ fainter around 
five-magnitude stars. Spectroscopy would be possible from flux 
ratios of 10"'*. Optimized interferometers with 10% transmission 
(Fig- El middle and bottom curves) would allow the detection of 
companions 10'* fainter around 5 to 7 magnitude stars and lO-' 
fainter around 8 to 10 magnitude stars. Spectroscopy would be 
possible from flux ratios of 10"^. 

Figure |4j) shows the baseline needed to resolve the primary 
star, that is, to reach the first visibility null for late dwarfs and 
giants. The horizontal line at 200 m corresponds to the maximum 
baseline of the VLTI, which imposes maximum H magnitudes of 
about 2.5 and 1.5 for MO dwarves and giants, respectively. Also 
the field of view of fiber-linked interferometers is limited to one 
Airy disk AjD, where D is the telescope diameter. Given that 
the baseline for reaching the first null is B = 0.6 1/i//?*, it gives 
a maximum recoverable separation of about BjD stellar radii, 
which is 100 stellar radii on the VLTI with 2 m telescopes. A 
baseline of 1000 m necessary for resolving a solar type star at a 
distance of lOpc would provide a field of 500 stellar radii. 

As specified in Sect. 13.2! the minimum error is attained in 
the photon noise regime when the working frequency intervals 
cover at least one period This implies, for the first null, a 
baseline variation AB of 



AB K 2- 



(12) 



At kilometric baselines, it represents a variation of 200 m for a 
distance of 10 stellar radii. The simplest way to achieve such 
variations is to use earth rotation with repositionable telescopes 
located at positions such that the baselines and orientations op- 
timize the observing time around the nulls. However, the double 
system may also be characterized at the price of a lower signal- 
to-noise ratio, even if the working frequency intervals cover 
much less than a period. Indeed, the companion of the spectro- 
scopic binary cr Puppis has been detected and characterized with 
only 3 snapshot observations of 5 minutes each; before, during 
and after the null, using the spectral dimension to increase the 
frequency coverage. 

4. Discussion 

4.1. Science cases 

We have demonstrated that a phase closure nulling experiment 
can bring astrophysically significant information on the close en- 
vironment of resolved stars. Figure |5] shows the lower limit of 
the number of targets available to phase closure nulling experi- 
ments as a function of basel ine length. The sam ple was obtained 
with the SearchCal utility ( Bonneau et al.ll2006l) and is limited to 
the stars whose apparent diameter can be reUably estimated from 
parallax and (spectro)photometric measurements available at the 
CDS. The minimum number of targets available for current ar- 
rays like the VLTI is ailOO in K and a!500 in J. The technique is 
not limited to observing unbalanced binary systems but can also 




Baseline Length (m) 



Fig. 5. Lower limit of the number of stars resolved at a given 
baseline in the J, H, and K bands. 



address a broader range of topics. These topics may change ei- 
ther by changing the nature of the resolved bright astrophysical 
target and/or the ofF-axis accompanying astrophysical signal. 

In star and planet formation, astronomers are interested in 
protoplanetary disks, which in general have the property of be- 
ing centro-symmetric so that the phase closure signal is zero. 
Any perturbating signal located off'-axis can then be detected by 
closure phase nulling: inhomogeneities on the disk surface, pres- 
ence of a forming planetesimal, presence of a collimated jet of 
outflowing wind partially screened by the disk. It also applies to 
other science cases holding disk-like cataclysmic variables, ac- 
tive galactic nuclei, proplanetary nebulae, gaseous disk around 
hot stars. The case of exoplanets is of particular interest, but the 
various kind of binaries must not be forgotten involving young, 
main sequence, or evolved stars and even the case of merging 
binaries. 

In a future where interferometry is sensitive to fainter targets, 
a possible application of this technique is the detection of the 
gravitational images due to a massive gravitational lens spatially 
resolved. In the following sections, we restrict ourselves to the 
companion of evolved and main sequence stars. 



4.1 .1 . Companions of evolved stars 

Clearly the simplest observational case, already permitted with 
current interferometers, is the study of close companions to 
giant stars. The scientific interest is twofold: improving our 
knowledge of the masses of giant stars and detecting plane- 
tary systems around earlier type stars than the ones known to 
date. It is indeed customary to look for companions of stars of 
masses > 1 .2R0, not in their main-sequence phase, where radial 
velocity measurements are rendered difficult by the paucity of 
spectral lines (furthermore broadened by the rapid stellar rota- 
tion), but in their giant, cooler phase. However, the detection of 
substellar companions around giant stars using radial velocities 
has always been problematic because the intrinsic pulsation of 
those stars would mimic the signature of a close substellar com- 
panion, and radial profiles can be altered by large, slowly rotat- 
ing, cool spots on the surface. Only a few detections have been 
reported to dat e ( Hatzes & CochranI Il998t iHatzes et al.1 120051: 
'Dollin ger et allllOOTT The phase closure nulling would provide 
a completely independent method of characterizing at least the 
closest companions, and give the orbit needed to measure indi- 
vidual masses. 
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We demonstrate this very possibility in a companion paper 
dPuvert et al .112009 ^. where we apply the phase closure nulling 
method for the detection of the close companion of the star 
Sigma Puppis and the characterization of this SBl binary sys- 
tem. 



4.1 .2. Companions of main-sequence stars 

These companions can range from brown dwarves (hereafter 
BDs) or hot Jupiters (HJs) to earth-like planets. With an expected 
flux ratio of a few 10"^ in the near infrared, brown dwarves could 
be detected with present-day telescope apertures at any distance 
from the star within the field-of-view permitted (xB/D). The ex- 
periment could bring two new pieces of information: the mass of 
the BDs (by solving the orbit) and the NIR spectrum, possibly 
different in close binaries than the spectra of isolated BDs. 

A few HJs, thus very close to their star (<20R*) and with 
temperatures of 1000-1500 K, will be detectable as well with 
a flux ra tio in the 10"*- 10'^ range in K band and longwards 
(see e.g. iFortnev et al.l l2008l) . A handful of HJs have already 
been photometrically detected with the IR AC and MIPS in- 
strument of the Spitzer S pace Telescope (ICharbonneau et al] 
12005 '; Knuts on et al. ]|2008') in transit experiments. Phase closure 
nulling experiments in this case would permit the NIR spectrum 
of the object, the orbit and thus the mass to be obtained, without 
being restricted to transit objects. Phase closure nulling exper- 
iments on dwarf stars with HJs would only be possible today 
on a handful of them since the baseline needed to null a Sun 
equivalent at lOpc is »1000m at 1.65jum. This restriction can 
be overcome in the near future as fibered i nterferometers wi th 
kilometer baselines are currently underway (iPerrin et al 

For any planetary companions still cooler and smaller than 
HJs, the flux ratio between the star and the companion rapidly 
becomes forbidding (>10'^) in the NI R and is m anageable only 
longwards of lOjum (see Fig. 2 in iTraub & Jucks 2002). In 
the N band, the above-mentioned constraint on the baseline 
length then becomes very challenging (>10'*m), and the signal 
transmission on such distances would be quite impossible for 
ground-based interferometers, unless techniques like the para- 
metric co nversion to vi sible or NIR wavelength, such as de- 
scribed in iDef Rio et al. 1 (120081) . are used before recombination. 



4.2. Towards a dedicated instrument 

We have discussed the main properties of phase closure nulling: 
around each null of the visibility function of the star, a spatial 
frequency interval exists where the signal of a companion be- 
comes stronger than any systematic error, hence potentially de- 
tectable. The errors on the companion detection, both in flux ra- 
tio and in position, depend roughly on the square root of the 
number of star photons detected, and the imageable region for 
a fiber-linked interferometer (the distance from the star where 
the companion's signature can be seen) is xB/D stellar radii. 
The need to be at a star null fixes a minimum baseline length. 
The need to reach a given sensitivity for the flux ratio implies 
either a given minimum size for the telescope apertures, thus a 
maximum imageable area, or/and a given integration time. The 
spectral resolution 'R does not need to be very high, typically 
xB/D. 

Today interferometers with baselines of a few hundred me- 
ters can already detect faint companions to the few closest giant 
stars using phase closure nulling. It seems possible to go further 
— and still use ground-based optical arrays — to characterize 



the faint companions, down to HJ masses, of the nearby main- 
sequence stars. This can be envisioned with an interferometer, 
that 

1. reaches the first zero of a solar-type star at lOpc in an in- 
frared band — already kilometric baselines in H band; 

2. has a spectral resolution of xB/D to enable detection of com- 
panions up to xB/D stellar radii. This is 500 for a 1000m 
baseline and 2m-class telescopes; 

3. permits long integration times on the zero crossing thanks to 
earth rotation and repositionable telescopes; 

4. has the highest attainable transmission: since our S/N only 
depends on the photon S/N of the primary, this should be as 
large as possible. Since we are only interested in phase clo- 
sure that is a robust measurement insensitive to relay optics 
wavefront disturbance or atmosphere stability, the attention 
should go on improving the transmission of the relay optics, 
not their optical quality. This is best realized with fully in- 
tegrated optics, fiber linking the foci of the telescopes at the 
primary focus. 

5. uses a few detectors with low noise for the de tection, such as 
a 4-pixel ABCD method. New developments dRothman et al.l 
12007) seem to imply that such low-noise IR detectors can be 
available in a few years from now. 

6. takes advantage of the new developments in integrated op- 
tics, permitting both beam recombination and spectrography, 
thus avoiding the losses at all the diopters of bulk optics 
used for beam transportation and in the spectrograph. Such 
a concept, der ived from on-going de velopments for visible 
spectrosco py (Le C oarer et al.l 1200 7*) is under investigation 
jKern et a l. 2009J and can be envision on the sky in the com- 
ing few years. 



5. Conclusion 

We have introduced both the concept of phase closure nulling 
and its application in the detection of faint close companions 
to resolved stars, accuratly measuring their position and spec- 
tra. This concept is the only one to date to permit such direct 
detection inside an Airy disk. We developed the complete error 
analysis of the phase closure nulling, and show how this method 
permits ground-based interferometry to complement and, in fa- 
vorable cases, compete with the direct detection methods con- 
sidered for future space-borne experiments. We discussed a few 
science cases that would benefit from phase closure nulling ex- 
periments and sketched the requirements for a new type of inter- 
ferometer dedicated to such studies. 



Appendix A: Error calculation 

A. 1 . Error on phase closure 

The formal error on the phase closure, in the presence of pho- 
ton noise and detector readout noise, was computed by IChellil 
( 1989) and Tatulli & Chelli (2005). Neglecting the coupHng be- 
tween the photon noise and the detector noise and retaining the 
dominant terms, the variance of the phase closure 0c may be ap- 
proximated by 



W2 

\ tel ^, 



1 



1 



1 



|V(«12)P 
tel" V 



2\ I 



|V(«23)P |V(«13)|2 
1 



|V(Ml2)F|V(M23)P|y(Mi3)P 
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1 



o-((f>c) dpj 



(with {pi,p2,P3} = {R*,r, s}). (A.2) 



The variances of the parameters are the diagonal elements of the 
matrix ['Ai x Af] ' . Let us first examamine the matrix 'Ai x Ai. 
It is given by 



1 50*30: 



1 50*50* . 

hk^-^ — hk ' ' 



1 .50* 2 



o-ldR^ ds 
cr} ds dr 



crjdR^ dr erf dr > 
* cr}dRi, ds * cr? ds dr * cr^^ ' 



ds 



(A.3) 



We neglect now the detector noise (photon noise regime), 
and we assume that we are working outside the null regions, 
so that the phase closure may be approximated by Eq. |7] Each 
element of the matrix above is the sum of nine terms. There are 
two types of terms 



^ /(Mh,)sin(27rMf„^)cos(2™*,^i) 

k 

Bin,pq = ^ g(Mf„)sin(2;rMf„s)sin(27rM*^j) 



(A.4) 
(A.5) 



where f{u) and g{u) are smoothly varying functions over the 
working frequency intervals Am/„. 

At last, we assume that the projected separation s is larger 
than a few stellar radii and that the working frequency intervals 
are larger than i"', but smaller than the mean frequencies: 



i > a few Ri, 



and 



M/„ s> Am/„ > s 



(A.6) 



Under these conditions, all the A terms are equal to zero and, 
given that ui2 + M23 + mi3 (non redundancy), the only non null 
B terms are those for which {/«) = \pq\- In addition, 

BinM = Yj ^(«L)sin^(2'rMf„s) * y X g(M/„) (A.7) 

k 

where ^(m/,,) is the average value of ^(m/,,) over the working fre- 
quency interval Am/„. The matrix 'AI x AI may then be written 



into the form 



where K is the total number of detected photoevents, A^tei the 
number of telescopes, Npi^ the number of pixels, and cr,. the de- 
tector readout noise. And |y(Mi2)|, I V(m23)|, |V(mi3)|, are the visi- 
biUty amplitudes of the observed object. 

A.2. Error on the parameters of a double system 

We want to estimate the errors on the parameters of a double sys- 
tem made of an extended uniform disk and a point source, from a 
least square fit of a set of Nc uncorrelated phase closure measure- 
ments obtained with an interferometer formed by A^tei telescopes. 
For simplicity, we assume that the telescopes are aligned, with 
a non redundant configuration, and that the direction of obser- 
vation varies little during the experiment and may be considered 
constant. Under these conditions, the object is characterized by 
3 parameters: the stellar angular radius the flux ratio r, and 

the projected separation s along the direction of observation. 

The Nc phase closures are sampled at frequencies Mf„{A: - Y(R-i,) 
I, Nc', In - 12, 23, 13), in frequency intervals Am/„. We define 
the Nc lines and 3 columns matrix Af, of element m jk, by 



'AI X AI ~ 

nL 



r'X(R^) rY(R^) 



rY{R^) 




(2nu,„axrYZ(R-^) 



(A.8) 



where K is the total number of photoevents detected during the 
whole observation and u„u,x the average value of the highest fre- 
quency. X{Ri,), Y(Ri,), and Z(Ri,) are given by 



1 



XiR*) 



dV^,(U|„) 



Vt(ui„)\ dR, 
1 



(A.9) 



Vl(u,„) 

1 dV^juln) 

Vl(u,„) dR^ 
1 



(A.IO) 



Vl(ui„) 



Z{R^) 



1 



(A. 11) 



Vl(u,„) 



The errors on the parameters of the system may now be com- 
puted. Inverting the matrix 'Ai x AI provides 



N^ 

2/n \ tel ^, 

o- (R*) ~ -J- X 
r^K 



1 



X(R,) - Y(R,) 



a-\r) 



(tHsIR^) 



N^ 
K 



X(R*) 



Ni 



X{R^) - YiR^) 
1 



(A.12) 
(A. 13) 
(A. 14) 



r^K (2nu,„axR*?Z{R^) ' 

Let us focus on the flux ratio and on the stellar separation. 

Numerical calculations show that the quantities X(R-i,)/[X{Ri,) - 
2 

Y(Ri,) ] and ZiR^,) are fairly constant and close to 1 in a wide 
range of experimental conditions, even for observations per- 
formed around the nulls. Then, in the photon noise regime, the 
errors on the flux ratio and the separation may be fairly approxi- 
mated by 



Me, 

o-(r) a — 



cr{s/R^) 



1 



Mel 



(A. 15) 
(A. 16) 
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